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In supersymmetric extensions of the particle physics Standard Model, gauge invariant combinations 
of squarks and sleptons (flat directions) can acquire large expectation values during a period of 
cosmological inflation. If the inflaton sector couples to matter fields via these flat directions, then 
new channels for efficient reheating, in particular via parametric resonance instabilities, are opened 
up. These can lead to efficient reheating induced by the flat directions even if the bare coupling 
constants are small. In this Letter we discuss various channels which yield this "enhanced reheating" 
effect, and we address some cosmological consequences. 

PACS numbers: 98.80.Cq, 11.30.Pb 
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I. INTRODUCTION 

There are strong indications that the early Universe 
might have undergone a phase of inflationary expan- 
sion [1] . The success of big bang cosmology puts severe 
constraints on early Universe cosmology, and therefore, 
if inflation is the right paradigm, then it must yield a 
way of understanding the origin of all observed matter, 
including the baryons of the Standard Model (SM). In 
spite of many attempts, it has proven hard to pin down 
the inflaton sector. Since its couplings with standard 
model fields in the bare Lagrangian are constrained to 
be very small, it is usually assumed that the inflaton is 
a SM gauge singlet, and in most models it is even con- 
sidered as an absolute gauge singlet (for reviews see e.g. 
[2,3]). This leads to the question of how the Universe 
reheats and produces the matter we observe today. 

It is commonly believed that the inflaton sector couples 
to either the SM sector or the matter sector of the min- 
imal extension of the SM, the Minimal Supersymmetric 
Standard Model (MSSM) through non-renormalizable in- 
teractions (for a review see [4]). In this context, however, 
it is not clear why the inflaton would decay at a sufficient 
rate into the SM degrees of freedom. This is the issue we 
will focus on in this Letter. We will concentrate on the 
case when the inflaton couples to either SM or MSSM 
fields via interactions which are suppressed by the scale 
at which new physics is being invoked, which could be 
either the GUT scale, string scale or the Planck scale. It 
turns out that the careful consideration of the evolution 
of the fields corresponding to the flat directions leads to 
a new mechanism of reheating which we call "enhanced 
reheating" 

As is well known, MSSM has flat directions, made up 
of gauge invariant combinations of squarks and sleptons 
(possibly including right-handed sneutrinos), which may 
acquire non- vanishing expectation values during inflation 



(see e.g. [2] for a review), thereby forming homogeneous 
condensates. The condensates may play a significant role 
in many cosmological phenomena (see e.g. [5] for recent 
reviews), such as generating the baryon asymmetry, or 
producing the dark matter particles by first fragmenting 
into Q-balls which then decay [6]. It has also been sug- 
gested that the cosmological density perturbations could 
be due to fluctuations of a MSSM flat direction conden- 
sate [7]. In this Letter we will study novel effects of 
MSSM flat direction condensates for reheating. 



II. ENHANCED PERTURB ATIVE REHEATING 

Even if at the level of the bare particle physics La- 
grangian there is no coupling between the inflaton and 
the MSSM matter fields, such a coupling will be gener- 
ated via gravitational effects. The resulting terms will 
typically correspond to a non-renormalizable field the- 
ory and will be suppressed by the Planck mass. We will 
model these terms via a simple superpotential term 



W D A— X*i* 2 , 
Mp 



(1) 



where <J> is the (singlet) inflaton superfield, X is the su- 
perfield containing the condensate, stands for MSSM 
matter superfields, Mp = 2.4 x 10 18 GeV is the reduced 
Planck mass, and A is a coefficient. Note that all what 
is needed is that X\&\& be a gauge invariant combina- 
tion of the MSSM fields. Then, regardless of any discrete 
or continuous global symmetry (e.g. R symmetry) that 
the theory may have, gravitational effects are expected 
to generate the above superpotential term. If the infla- 
ton is a gauge singlet up to Mp, A ~ 0(1) typically. On 
the other hand, if $ is non-singlet under new physics at 
a scale M now < Mp, A will generally be suppressed by 
powers of M new /Mp. 
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Around the minimum of its potential, the inflaton su- 
perpotential can be approximated by m^ 2 /2. At the 
end of inflation (X) = Xi and, treating the inflaton as 
a real scalar field, eq. (f ) results in a Lagrangian of the 
form 



(2) 



where tp is a Dirac fermion containing the fermionic com- 
ponents of and 4 , 2 ) and \ denotes the scalar compo- 
nents of the ^ superfields. 

A phase of primordial inflation driven by the slowly 
rolling scalar field </> sweeps the Universe clean of any 
particle excitations which might have been present be- 
fore inflation, but leads to vacuum fluctuations which 
generically will produce quasi-homogeneous expectation 
values (condensates) of light fields such as Xj. After 
the slow-roll approximation breaks down, the inflaton 
field will begin a phase of coherent oscillations about 
the minimum of its potential. During this phase, the 
coherent oscillations will dominate the energy density of 
the Universe. Due to couplings of the inflaton field to 
other fields, the oscillations will transfer some of their 
energy density to the matter particles. Upon thcrmaliza- 
tion of these particle excitations, a radiation-dominated 
Friedmann-Robertson- Walker (FRW) universe will re- 
sult, as in the hot big-bang model. This process is 
called reheating. Successful Big Bang Nucleosynthesis 
(BBN) [8] requires that the reheat temperature of the 
universe obey T R > O(MeV). 

In the initial analyses [9] , the energy transfer from the 
coherent inflaton field to matter was treated perturba- 
tively, using the lowest order analysis to determine the 
rate Td with which a single cf) particle at rest decays into 
quanta of the fields which <fi couple to. While the Hubble 
constant H exceeds the decay rate, the inflaton will lose 
energy primarily to the expansion of space, but as soon 
as H drops Td, the energy of cf) goes primarily into matter 
particles, resulting in a gas of matter with a temperature 
T R given by 



o.i(r d M P ) 1/2 . 



(3) 



In our toy model (2), the inflaton can decay perturba- 
tively via several channels. In the limit that the inflaton 
mass ttt.0 is much larger than the mass of ip and \, the 
perturbative decay rate of (j> into fermion and scalar pairs 
is given by 



Tt = A 



2 m <t> 
47T 



(4) 



Let us assume that the perturbative decay rate in the ab- 
sence of the X condensate is T (with associated reheat 
temperature Tr ,q determined via (3)). From (4) it imme- 
diately follows that a large value of Xi leads to enhanced 



perturbative reheating. The condition for enhanced re- 
heating (i.e. Ti > r ) is 



A 



> 10 3 



T 2 
m&Mp 



(5) 



As an example, consider the chaotic inflation model with 
to = 10 13 GcV and T R . = I0 9 GcV. Then the bound 
in (5) translates to AAi > 10 13 GeV. 



III. ENHANCED PARAMETRIC REHEATING 

As first discussed in [10] (see also [11]), in many mod- 
els there is a non-perturbative and much more efficient 
reheating mechanism which makes use of the parametric 
resonance instability in the field equation of the particles 
being produced (ip and x m our model) in the presence 
of an oscillating inflaton field. Note that, as first dis- 
cussed heuristically in [10] and then analytically in [12] 
and [13], the instability is present even when taking the 
expansion of space and the resulting decay of the infla- 
ton oscillation amplitude into account. This mechanism 
leads to an initially highly non-thermal distribution of 
particles (including inflatons), while full thermal equilib- 
rium takes much longer to establish (see e.g. [14] for a 
recent study). 

The Lagrangian (2) contains three separate channels 
for resonance. The last term on the right-hand side of (2) 
corresponds to two-particle induced resonant production 
of x bosons and was treated in detail in [15]. The second 
term leads to single-particle induced resonant produc- 
tion of x bosons and was analyzed in [13], while the first 
term lead to parametric production of ip fermions and 
was treated in detail in [16]. The third term dominates 
for large values of Xj. 

Let us first focus on the third term on the right hand 
side of (2). If we neglect the other two terms, then in the 
presence of the oscillating inflaton field, the equation of 
motion for x will lead to exponentially growing solutions 
with a Floquet exponent /i much larger than the Hubble 
expansion rate. The instability occurs for a large range 
of comoving momenta (and thus the instability is called 
of broad resonance type [15]). The instability is effective 
as long as [15] 



A( 



(6) 



The effective coupling constant during inflation is 
bounded by the amplitude of the induced CMB 
anisotropies, leading to the constraint 



A 2 (^) 2 < HT 6 . 
y Mp> ~ 



(7) 



making use of the fact that an interaction coupling con- 
stant arising in (2) will induce a one-loop correction to 
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the quartic inflaton self coupling which is bounded by the 
standard fluctuation analysis (see e.g. [17] for a review). 

During inflation, long wavelength quantum fluctua- 
tions of X accumulate in a coherent state with a maxi- 
mum quasi-homogeneous field value given by [3] (X 2 ) = 
3Hf /8TT 2 m 2 K , which for condensate masses smaller than 
the Hubble constant during inflation can be many orders 
of magnitude larger than Hi. The exception is when su- 
persymmetry breaking by the inflaton energy density [18] 
results in a correction « +H 2 to the (mass) 2 of X dur- 
ing inflation.* Another way to obtain a large expectation 
value Xj is if X is exactly massless during inflation and 
obtains a non-vanishing mass during a phase transition 
which takes place during the early stages of reheating, in 
analogy to how the QCD axion field develops an expec- 
tation value after the QCD phase transition. Note that 
for field values larger than Mp, gravitational corrections 
to the potential may not be under control, and hence we 
will consider Hi <C Xi < Mp at the end of inflation. 

The scenario of enhanced parametric resonance is now 
as follows. At the end of the slow-rolling phase of in- 
flation, <j> ~ M P , and, since the inflaton mass is 
constrained by the amplitude of the CMB fluctuations 
to satisfy < 10 _6 Mp, the broad resonance condition 
(6) is satisfied and resonant \ production occurs. As 
a consequence of the energy loss to particle production 
and of the expansion of space, the amplitude of </> oscilla- 
tions will decrease. If the condition (6) is satisfied when 
resonant particle production becomes explosive, inflaton 
oscillations will decay in the broad resonance regime [15]. 
When the condition (6) ceases to be satisfied the broad 
resonant decay shuts off. 

Even after the period of broad resonance ends, para- 
metric resonance via the first two terms in (2) will con- 
tinue for a while. To see this, focus on the resonant 
production of x particles via the single <j> coupling in (2). 
The perturbative decay rate is the same as the perturba- 
tive decay rate (4) into fermions via the first term on the 
right hand side of (2). As studied in detail in [13], the 
parametric resonant decay of </> proceeds by exciting par- 
ticles in narrow resonance bands, the first at k = m^/2. 
Decay into the lowest resonance band dominates the over- 
all decay rate, and, focusing on this band, the decay rate 
is given by 



Hm 2 
16ir 2 4>'- 



:Sinh *[X 2 (fLf 



M P ' Hm, 



(8) 



(see Eq. (125) of [13]). 

If the argument of the hyperbolic function is smaller 
than 1, the result reduces to the perturbative decay rate. 
Thus, this channel has enhanced parametric reheating 



*Note that a correction « — H 2 typically results in Xj 3> Hj, 
since X will in this case settle at the minimum of the potential 
which is far away from the origin. 



provided the argument is larger than 1 , which is the case 
if (taking into account the fact that the inflaton potential 
is taken to be quadratic) 



A" 



Mpmj 



> 1. 



(9) 



By inserting the saturation values of (6) and (7) into (9), 
it follows immediately that this condition remains sat- 
isfied once the broad resonant decay shuts off, provided 
that M P > lO 3 m0. 

Once the value of H drops below mx, the X con- 
densate will also begin damped oscillations. Thus, the 
conditions in (6), respectively (9), may cease to be sat- 
isfied before broad, respectively narrow, resonant decay 
completes. However, the oscillations of X lead to an- 
other channel of parametric reheating, namely reheating 
induced by the coupling of X to other matter fields in 
the part of the Lagrangian which is independent of the 
inflaton (see [19] for a discussion of this scenario). 

To conclude this section, we have demonstrated that 
large expectation values of scalar condensates such as 
arise in the MSSM open up several channels for enhanced 
reheating via parametric resonance. 



IV. SOME CONSEQUENCES 

Let us first discuss the consequences of the above anal- 
ysis for the "temperature" Tr after reheating, by which 
we mean the fourth root of the energy density in matter 
after the inflaton has lost a majority of its energy into 
matter t There are various possible scenarios. 

In the first case, the broad resonance condition (6) is 
satisfied when the inflaton begins to oscillate, which oc- 
curs when <f> ~ Mp. In this case, the energy transfer into 
matter happens explosively via broad parametric reso- 
nance, yielding the value (remembering again that the 
potential for <p is quadratic) 



(m^Mp 



,V2 



(10) 



which, for values ~ lCP 6 Mp obtained from the am- 
plitude of the observed CMB spectrum under the hypoth- 
esis that adiabatic perturbations of the inflaton generate 
the observed anisotropies, is about 10 16 GeV. This sce- 
nario is realized quite naturally if the initial value of the 
condensate is close to the Planck scale. 



Hf the state after this stage of "preheating" [12] were ther- 
mal, then the quantity Tr thus computed would indeed be 
the usual thermodynamical temperature of matter. However, 
if the reheating proceeds via parametric resonance, then the 
state of matter immediately after preheating will be far from 
thermal. 
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If the broad resonance condition (6) is not satisfied 
when the inflaton starts oscillating, or ceases to be satis- 
fied before broad resonant decay completes, but the con- 
dition (9) for effectiveness of the narrow resonance insta- 
bility via the second interaction term on the right hand 
side of (2) is satisfied, then once again the energy trans- 
fer to matter will be immediate and the "temperature" 
after the preheating will be given by (10). 

The third scenario arises when neither (6) nor (9) are 
satisfied, but when the perturbative decay rate in the 
presence of the condensate is larger than m x , the 
mass of the condensate. In this case, reheating completes 
before X starts its oscillations. In this case, the reheat 
temperature is given by (3), with Td = Fr, yielding 

T R . 3 x 10-^(11) V2 . (11) 

Finally, if none of the above resonance conditions are 
satisfied and in addition T\ < mx, then the condensate 
will begin to oscillate once the Hubble parameter H drops 
below mx , and this will lead to new production channels 
for the production of matter via X decay. Assuming that 
the decay is rapid on a Hubble time scale, the resulting re- 
heat "temperature" Tr (quotation makes indicating once 
again that the quantity is simply a measure of the energy 
density of the matter particles after the decay of the X 
condensate) will be given by 

T R ^ (m x Xj) 1/2 , (12) 

unless the result is smaller than what is obtained from 
(11), in which case most of the energy of matter particles 
results from the inflaton decay which starts once T\ = H , 
i.e. after the condensate has started to oscillate. 

We have seen that reheating can be efficient even for 
very small inflaton couplings to matter. This can have 
good or bad effects for early Universe cosmology. On 
the positive side, the enhancement of reheating by the 
X condensate may make it possible to obtain reheating 
temperatures sufficiently high for successful baryogenesis 
and/or dark matter production even when the reheat- 
ing temperature Tr j0 (calculated using the perturbative 
decay rate in the absence of the condensate) is too low. 
However, enhanced reheating should not lead to the grav- 
itino problem [20]. If T R > 10 9 GeV, then there will be 
a gravitino overproduction problem unless the gravitino 
mass TB3/2 obeys 1^3/2 > 20 TeV (this happens, for ex- 
ample, in models of anomaly-mediated supersymmctry 
breaking [21]). Otherwise, late entropy generation will 
be required to dilute the excess of gravitinos. As we have 
seen above, if the parametric resonance channels for re- 
heating are open as a consequence of a large expecta- 
tion value of X, then Tr can be of the order 10 16 GeV 
and we will be faced with the gravitino problem. Even if 
the parametric resonance channels are closed, then, using 
(11), we conclude that there will be a gravitino problem 
if 



AX/ > 10~ 5 M P . (13) 

A second stage of inflaton like in thermal inflation [22] 
can dilute the gravitino density and thus solve the grav- 
itino problem. If Xj > Mp at the end of inflation, then 
such a second stage of inflation might naturally arise in 
our framework: once H drops below the value of mx 
and the condensate can start to move, it will initially 
be slowly rolling and thus generate a second period of 
inflation. 

Another consequence of the enhanced reheating sce- 
nario discussed here concerns the origin of cosmological 
fluctuations. As is well known quantum fluctuations of 
the X condensate will lead to isocurvature fluctuations 
as in the case of axions (see e.g. [23]) which upon the 
decay of the condensate will convert to adiabatic fluctu- 
ations via the curvaton mechanism [24]. However, the 
X fluctuations will also lead to a spatial modulation of 
the reheating rate on super-Hubble scales and will thus 
generate fluctuations via the modulated decay mechanism 
[25] (this has already been pointed out in the context of 
the MSSM in [26]. As discussed in [27], the fluctuations 
which are generated in this way are primordial isocur- 
vature fluctuations, and their production is completely 
consistent with causality. 

The situation will be different when 

r < m x , (14) 

Xj < Mp , and if the parametric resonance channels for 
enhanced reheating are not open. In this case, the X 
field starts oscillating at H = mx , before reheating com- 
pletes. From then on, (X) oc H leading to Tj(t) oc H 2 . 
This implies that reheating will remain inefficient until 
Td settles at its value in the minimum Tq. The reheat 
temperature of the universe will then be Tufi- 

An important issue to take into consideration is early 
oscillations of X due to thermal effects from perturbative 
reheating [28]. The X field has gauge and/or Yukawa 
couplings (collectively denoted by h) to other fields. If 
hX\ < T, these degrees of freedom will be in equilib- 
rium with the instantaneous thermal bath from reheat- 
ing. This in turn induces a thermal mass hT for the X 
field which will trigger its early oscillations if hT > H . 
Early oscillations, if they start before reheating com- 
pletes, would render the reheating inefficient as noted 
above. There will be no early oscillations if hXj > T or 
hT < H. This requires sufficiently large or small values of 
h, respectively, which can naturally occur for the MSSM 
flat directions and right-handed sneutrinos respectively. 
The early oscillations will not be important in the case 
of parametric reheating. In this case, fields which are 
coupled to X via gauge and/or Yukawa couplings reach 
thermal equilibrium on time scales much longer than that 
of the resonant inflaton decay [14]. 
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V. CONCLUSIONS 



project. AM thanks the ASPEN center for physics for 
their kind hospitality during the course of this work. 



In this Letter we have discussed enhanced reheating 
channels which arise as a consequence of the presence of 
scalar condensate fields which acquire a large expecta- 
tion value Xj during inflation. The scalar condensates 
can result in large effective couplings of the inflaton to 
matter, which can enhance both the perturbative decay 
rate and the efficiency of parametric resonance instabil- 
ities. We have seen that efficient reheating is therefore 
possible even if inflaton has only gravitational couplings 
to matter fields in the bare Lagrangian. 

Enhanced reheating can have many consequences for 
cosmology. Generically, the reheating temperature will 
be higher than what would result in the absence of the 
condensates. In supersymmetric models this can make 
the gravitino problem worse. 

Since scalar condensates acquire fluctuations during in- 
flation which are super-Hubble-scale at the end of infla- 
tion, they provide a realization of the modulated fluctu- 
ation scenario of [25]: at the time of reheating, super- 
Hubble scale entropy perturbations are generated as a 
consequence of the space-dependent effective coupling 
constants. This will relax the bounds from density per- 
turbations on the inflaton mass m<f,. 

In summary, scalar condensates can lead to a success- 
ful scenario of chaotic inflation where the inflaton has 
only gravitational couplings and <C 10 13 GeV. Im- 
mediately after the phase of slow-roll inflation comes to 
an end and the inflaton field begins to oscillate, the large 
amplitude of the condensate Xi leads to large effective 
couplings in the dimension 5 and 6 operators which cou- 
ple the inflaton, the condensate and matter fields. This 
opens up channels for decay of the inflaton via broad 
parametric resonance which would be closed in the ab- 
sence of X. Thus, rapid reheating takes place. The de- 
cay will produce massive particles which can provide the 
dark matter and/or (in the case of right-handed neutri- 
nos) the lepton asymmetry required for leptogenesis. At 
the same time, the super-Hubble-scale fluctuations of the 
condensate lead to isocurvature fluctuations both of ax- 
ionic type and of the type generated through modulated 
decay. A detailed study of enhanced reheating and its 
consequences for baryogenesis and dark matter will ap- 
pear in a separate publication. 
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